Reversal of Nonlocal Vortex Motion in the Regime of Strong Nonequilibrium 
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We investigate nonlocal vortex motion in weakly pinning a-NbGe nanostructures, which is driven 
by a transport current / and remotely detected as a nonlocal voltage V n i- At high /, the measured 
V n i exhibits dramatic sign reversals that at low and high temperatures T occur for opposite polarities 
of J. The sign of V n i becomes independent of that of the drive current at large \I\. These unusual 
effects can be nearly quantitatively explained by a novel enhancement of magnetization, arising 
from a nonequilibrium distribution of quasiparticles at high T, and a Nernst-like effect resulting 
from local electron heating at low T. 

PACS numbers: 74.25.Qt,74.25.Fy,74.78.Db,74.78.Na 
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Motion of the Abrikosov vortex lattice in type-II su- 
perconductors results in strong deviations of the quasi- 
particle distribution function from that in equilibrium 

esq 

when the lattice is strongly driven by a trans- 
port current. Close to the critical temperature T c , over- 
heating of quasiparticles within the vortex cores leads to 
a shrinkage of the cores, accompanied by decreasing the 
effective viscosity coefficient r\ - the Larkin-Ovchimiikov 
(LO) instability 0], while the quasiparticles outside the 
cores remain in thermal equilibrium. At low T, the entire 
quasiparticle subsystem is heated because of the larger 
electron-phonon collision time. This results in an expan- 
sion of the cores instead of their shrinkage, while r\ again 
decreases In both cases, the current-voltage [V^(/)] 

characteristics are very nonlinear - can become even hys- 
teretic Q - and are in fact so similar that the difference 
can be resolved only via a quantitative analysis 0, 
However, vortex shrinkage and vortex expansion are dif- 
ferent effects and should lead to qualitative differences in 
other properties. 

In this Letter, we report novel effect in the recently 
discovered nonlocal vortex flow in the transversal flux 
transformer geometry (TFTE) 0, 0] , which allow a clear 
distinction of the above two opposite types of nonequilib- 
rium. We apply a drive current / in one part of the sam- 
ple (local lead) and measure the voltage response (V n i) 
in a remote part of the superconductor connected with 
first one via a channel of the same material [see the inset 
to Fig. [Ua)] ■ I n such a geometry, one can probe changes 
in the vortex lattice which occur in the local lead via 
changes in the interaction between vortices in the local 
lead and vortices in the rest of the sample. In this way, 
we can detect a novel nonequilibrium enhancement of the 
magnetization of the superconductor in the LO state with 
respect to the equilibrium magnetization and observe a 
Nernst-like signal at low T. 

Previously, V n i(I) was investigated in the linear re- 



sponse regime [H, The main features of these stud- 
ies can be accounted for by a simple model of locally 
driven vortices pressurizing those in the channel by re- 
pulsive vortex- vortex interaction / applied between 
the contacts 1 and 2 in the inset to Fig. [TJa) decreases 
exponentially in the perpendicular channel, with a decay 
length W/ir < L [& if. Thus, « n<f> WX driven vor- 
tices face w n^WL vortices in the channel, where X is 
the effective length over which the driving force fd r (per 
unit vortex length d) acts, = B/<j)Q the vortex den- 
sity, </>o the magnetic flux quantum, B = B ext + imjM, 
B ext the external magnetic field, M the magnetization, 
and hq — An • 10 _7 Vs/Am. The driven vortices push 
or pull those in the channel by exerting a pressure p = 
{n c f > WX){fdr/W). The resulting force pWd is balanced 
by the total frictional force (n^W L)(rju n id) on the vor- 
tices in the channel (which move at velocity u n {). For a 
superconductor with a large magnetic penetration depth 
A, i.e., 77,0 w B ext /<p , using V nl = WB ext u n i for the 
voltage detected at the probes 3 and 4, one obtains 



Vm = WB ext Xf dr / V L 



(1) 



At low /, i.e., close to equilibrium, fd r is given by 
the Lorentz force fi = j<hu where j is the transport 
current density. In Ref. Q, X = W led to V n i = 
(WB ext cj) /r]Ld)I = R n il. This reproduced the observed 
V n i cx / and V n \ oc 1/L even in the presence of pinning 

m 

Our d = 40 nm thick a-Nbo.7Geo.3 samples are pro- 
duced by electron-beam lithography and magnetron sput- 
tering. The local current leads (1,2) are connected to the 
nonlocal voltage probes (3,4) via a perpendicular channel 
of L = 2 /Ltm and W = 250 nm. All data for Vi(I) refer 
to passing / between 1 and 3, and measuring Vi between 
2 and 4. Since W is also the width of all other narrow 
sample parts, in particular that linking 1 and 2, Vi(I) 
and V n i(I) can be compared directly. Measurements of 
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FIG. 1: Typical local (black lines) and nonlocal (colored sym- 
bols) V(I) curves, (a) T = 0.75 K = 0.26 T c (at B ext = 3.0 T, 
6 = 0.64), and (b) T = 2.50 K = 0.85 T c (at B ext = 0.50 T, 
6 = 0.50). Inset to (a): Sample geometry; L — 2 /zm, 
W = 250 nm. Inset to (b): Saturation voltages of V n i plotted 
against b, for low (red) and high (blue) T. 



V(I) provided all relevant parameters of our samples: 
T c = 2.94 K, the normal-state resistivity p n = 1.82 /j.Slm, 
-{dB c2 /dT) T=Tc = 2.3 T/K, where B c2 is the equilib- 
rium upper critical magnetic field, and the Ginzburg- 
Landau (GL) parameters K = 72, £(0) = 7.0 nm, and 
A(0) = 825 nm. The low pinning in a-Nbo.7Geo.3 allowed 
for dc measurements of V n i ~ 10 — 200 nV, which was at 
the level of R n i ~ 0.1 ft in the low-/ linear regime. All 
measurements were carried out in a 3 Hc cryostat, with 
Bext perpendicular to the film plane. 

Typical results for the two limiting cases of low (T = 
0.75 K = 0.26 T c ) and high (T = 2.50 K = 0.85 T c ) tem- 
peratures are shown in Fig. [lja) and Fig. [IJb), respec- 
tively. The Vi(I) curves exhibit a nonlinear shape char- 
acteristic of strong-nonequilibrium (SNEQ), originating 
either in (a) electron heating 0, 0, [f| or (b) LO vortex- 
core shrinking [3, 0, Si- On the other hand, V n i(I) dis- 
plays the previously observed linear, antisymmetric de- 
pendence [i.e., V n i(— I) = —Vni (I)] only at low I. Upon 
increasing /, sudden sign reversals of V n i are observed 
in both regimes: at a certain J, the antisymmetric sig- 
nal converts into a symmetric one. The sign of V n i can 
be unambiguously attributed to the following directions 
in the inset to Fig. [IJa): at low positive (negative) /, 
the positive (negative) V n i corresponds to vortex motion 
upwards (downwards) in the channel. When / is high, 



vortices move either downwards (T <C T c , V n i < 0), or 
upwards (T — > T c , V n i > 0), irrespective of the direction 
of I. The saturation values of \V n i\ at high I are plot- 
ted vs b — B ext /B c2 in the inset to Fig. [ljb). In both 
cases, nonzero values are observed only at intermediate 
b, with a maximum efficiency around b = 0.6 (6 = 0.45) 
at low (high) T, similarly to thepreviously observed B^r t 
sweep traces of V n i at low I d, 0| ■ As argued in Ref . 0] , 
the vanishing of V n i at low B ext is presumably related to 
becoming smaller than the density of pinning sites, 
whereas V n i(B ext — > B c2 ) — > because the sample goes 
to the normal state. 

We first discuss the regime T -C T c . Assigning the 
corresponding high-j SNEQ state to electron heating to 
T = T* above the bath temperature T was successful 
in explaining the measured VAI) of Refs. @, 0, 0|- An 
analysis of the present Vi (I) [8j within the same frame- 
work permits to extract T*(Vj) and, using V(I), also 
T*(J), which is more convenient for a comparison with 
the V n i(I) data (see below). The hot electrons penetrate 
into the channel, which remains at T = To, roughly up to 
L T = \fDT\) ~ 295 nm - W. Here, D = 4.80 • 10~ 5 m 2 /s 
is the diffusion constant, and to ~ 1.82 ns is the relax- 
ation time of the hot electrons, resulting from the men- 
tioned analysis Q. Hence, there is a T gradient which 
leads to a thermal driving force f-r = — S^VT and conse- 
quently to the Nernst effect. is the vortex transport 
entropy Q. The Nernst effect should lead to vortex mo- 
tion downwards, which agrees with the observed V n i < 0. 
Since T* - Tq ~ 1 K typically, the observed temperature 
gradients |VT| - (T* - T )/L T ~ 1 K/fan are much 
larger than in usual measurements of the Nernst effect. 

The above is elaborated in Fig. [21 where the result 
for B ext = 3.0 T (b = 0.64) is analyzed more closely 
The shape of V n i{I) in Fig.QJa) suggests to consider the 
symmetric (+) and antisymmetric (— ) parts of V n i sepa- 
rately via V^j(I) = [V n i(I) ± V n i(—I)]/2, which is shown 
in Fig. EJJa). V~ t (I) at low / is fairly linear as expected, 
since fd r = /l, while V^(I) is very small. Upon increas- 
ing /, this is followed by a rapid suppression of V~, (I) and 
a simultaneous growth of V^ (I) < to a constant value 
comparable to that of the maximum V~i(I) > 0. Return- 
ing to Fig. IDJa), one can note that this dramatic change 
occurs around I where V(I) ~ R n I, signifying the tran- 
sition to the normal state in the local region 0, 0] and 
consequent vanishing of f^. Furthermore, |7| where the 
sign of V n i changes steeply on the I > side (ft and 
act oppositely) coincides with |/| where V n i has a local 
minimum on the I < side (Jl and fr add); in both 
cases, this marks that only fr remains effective at higher 

In the main panel of Fig. [2jb) , we plot T* (I) extracted 
according to the electron heating model [1, 0| in the su- 
perconductingstate and from noise measurements in the 
normal state [8[, whereas in the inset we show a sketch 
of the T profile along the sample. One can see that the 
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FIG. 2: (a) Measured V n i at T <C T c : antisymmetric (red) and 
symmetric (black) part of the nonlocal signal at T — 0.75 K 
and Bext = 3.0 T. (b) Effective electron temperature T*(7), 
where the black line stems from an analysis of Vi(I), and the 
red line from noise measurements in the normal state. Inset: 
Sketch of the temperature profile along the channel. 



electron heating is basically absent at low 7, then sets 
in very steeply until it reaches T c (B ext ) that represents 
B C 2 {T) 0, [f| Q , after which it changes with 7 only weakly. 
The nearly flat V^(I) at high 7 hence corresponds to 
T* » T c {B ext ), so \VT\ » [T c (B ext ) - T Q }/L T = ST/L T . 
Using Eq. |T]), we can extract from our data by focus- 
ing on the saturating values of Vjj (7 ) ■ We approximate 
fdr = /t ~ S^ST/Lt and X f=a Lt to obtain 



Sj, = V n i4>o/ R n i5Td 



(2) 



which does not contain Lt- Since S</, and R n i depend 
on the properties of the channel (where T = To), the 
observed Kit (I) sw const, follows straightforwardly. In 
the (B ext ,T) range of our data, we find Sj, ~ 0.1 — 1.5 • 
10~ 12 Jm _1 K _1 [8[, which is in reasonable agreement 
with a theoretical estimate ~ 0.1 — 0.2 • 10~ 12 Jm _1 K _1 
obtained by using the Maki formula [ill ], as well as 
with experimental data on films of Nb (0.05 — 1.5 • 10~ 12 
Jm-iR" 1 ) [13 and of Pb-In (0.2 - 5 ■ 10~ 12 Jm"^" 1 ) 

m ~ n 

We now turn to the regime T — > T c . An analysis [8j] 
of the Vi(I) in the spirit of Refs. @, d reveals that this 
SNEQ state corresponds to the LO vortex-core shrink- 
ing [l| , with TkTo everywhere because electron heating 
is strongly suppressed close to T c 0, 0, @. Vi(I) for 
T <C T c and T — > T c are at first glance rather similar, so 
the difference becomes obvious only through a numeri- 
cal analysis 0, In contrast, the qualitatively different 
V n i(I) curves in Fig. 1 leave no doubt that we are deal- 
ing with two distinct SNEQ phenomena. As before, the 
shape of V n i(I) [see Fig.QJb)] suggests to consider V^(I) 
and V n 7(/) separately, which is shown in Figs. EUa) and 
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FIG. 3: Measured V^(7) (a), and V~(I) (b) at T = 2.50 K = 
0.85 T c and 6 = 0.45 (black), 0.50 (red), 0.55 (green) and 0.60 
(blue), (c) Calculated j s /jaL and |A|/Ao (inset) vs r/£(u = 
0) for different u/u LO - (d) Calculated (M neq - M eq )/M eq 
against u/ulo at different b (as indicated). 



|3Jb), respectively. V^7(7) at low 7 is linear for small b, 
which implies the presence of /z, , whereas this is difficult 
to claim for higher b where the signal is small over the 
entire 7 range. At high 7, however, V~i(I) is small re- 
gardless of b. V^"(J), on the other hand, increases with 
increasing 7, and eventually saturates at a value compa- 
rable to that of the maximum Nernst signal at low T, 
albeit with the opposite sign. The smallness of V^J(J) 
at high 7 suggests inefficiency of in this regime. This 
can be understood by recalling [see Fig. [ljb)] that the 
Vi (I) for these 7 is close to the normal-state dissipation, 
which means that most of the current is normal |l| - and 
normal current does not contribute to ]l- 

Since /x is negligible and T w To, there must be 
yet another driving force which governs the TFTE at 
high 7. Below we show that this force has the same 
origin as the LO effect on V/(7), that is, a deviation 
Sg(e) of the quasiparticle distribution function <?(e) from 
ffeg(e) = tanh(e/2fcsT) = g(e) —Sg(e) in equilibrium. An 
additional consequence of Sg is an enhancement of the 
supercurrent density j s flowing around the vortex core, 
which can be calculated following P, [3] 



7 |A|5 ff (|A|) 



2e 



p n e \4A:b7 c 

'(3) 

where A = |A|exp(i<y9) is the order parameter and A 
the vector potential. The term oc |A| 2 corresponds to 
the equilibrium contribution to j s in the GL model, and 
the term oc Sg to the SNEQ correction. Sg is positive 
for energies less than the maximal value |A| TO(M of the 
order parameter in a single- vortex cell [H, and |A| is en- 
hanced near the vortex core [see the inset to Fig. [3jc)]. 
Both these factors lead to a growth of j s near the vor- 
tex core [see Eq. (3)]. Therefore, the magnetic moment 
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m = (1/2) J[r x j s ]dS ce ii of each cell in the vortex lattice 
increases in the LO state. 

We base our model on addressing m to find M along 
the direction of B ext , which is alternative (but much sim- 
pler with regard to the role of Sg) to using the Gibbs free 
energy density for the same purpose. Qualitatively, m 
of a given cell creates a dipole magnetic field which in 
the surrounding cells opposes B ext , hence an increase of 
j s results in a stronger diamagnetic response. Note that 
the same argument can be used to explain increase of 
the equilibrium diamagnetism of the mixed state as T 
decreases. Quantitatively, we have to determine 17(e) and 
|A|. We follow the LO model and solve numerically the 
modified GL equation for |A| (sec Eq. (A49) in [Uncou- 
pled with the equation for g(e) (see Eq. (A45) in lj). 

In Fig. EJc), we plot exemplary j s /jcL vs re- 
duced radial coordinate r/£, where jgl — 0.93A (1 — 
T '/Tc) 1 / 2 /£p„e is the equilibrium GL depairing current 
density, A ~ 3.06fc B T c (l - T/T^ 1 / 2 , and £ corresponds 
to that at zero vortex velocity u. Results are shown for 
three different u relative to the LO vortex velocity ulo 
[H; the corresponding |A|/Ao is shown in the inset by 
the same colors. By summing up the resulting m of each 
cell, one can find the difference 8M = M neq — M eq of the 
noncquilibrium (M neq ) and equilibrium (M eq ) magneti- 
zation. This is presented in Fig. Efd) . The maximum of 
6M/M eq occurs for b ~ 0.2 (at u/ur,o ~ 1). At smaller 
b, the enhancement of j s near the core gives a small con- 
tribution to m. At larger b, the suppression of |A| at the 
cell boundary [due to Sg(e) < for e > |A| max ] becomes 
important. We show results up to u = 2ulo, where the 
LO approach becomes invalid at T ~ 0.85 T c . 

The spatial variation of M across the boundary be- 
tween the local region and the channel occurs over a 
length of about the intervortex distance do ~ \J 4>o/ B ex t, 
and induces a current density jm = V x M that flows 
along that boundary. This current creates a force /m = 
Jm^o that is again independent of the direction of /, pulls 
the vortices toward the local lead (which results in V n i > 
0), and dominates the total f^ r in the SNEQ regime near 
T c . The typical \6M\ ~ \M eq \ ~ {B c2 - B ext )/2p K 2 ~ 
35 A/m (= 88 ^T at B ext = 0.45 T) is rather small but 
appears over a very small distance ao(B ext = 0.45 T) ~ 
70 nm, thus providing jm — 500 M A/m 2 which is of the 
same order as the transport current densities we used - 
as / = 1 //A corresponds to j = 100 M A/m 2 . We again 
employ Eq. (1) to estimate V n i- Since j'm = dM/dx « 
\5M\/a$ and X w ao, with f c i r = f^j we obtain 



V nl = [WB ext a /r]L}j M cj) = R nl \5M\d 



(4) 



from which ao has dropped out again. Inserting typical 
values of R n i w 0.1 fl and |<5M| « 35 A/m, we find V n i R* 
140 nV, which is quite close to the measured values. 

In view of the simplicity of our model, the agreement 
between the experiment and theory is rather remarkable. 



A full quantitative account of the phenomenon would re- 
quire inclusion of other effects on the interface of the local 
region and the channel - such as details of entry/exit tra- 
jectories for the fast and slow vortices, etc. However, 
these corrections may depend on the sample geometry, 
and we believe that the main physics of TFTE close to 
T c is captured by our model. 

In conclusion, nonlocal measurements allowed us to 
qualitatively distinguish two different types of vortex mo- 
tion in strong nonequilibrium. According to our theory, 
close to T c a new type of nonequilibrium magnetization is 
built up in the drive wire, which pulls the vortices towards 
the drive channel. At low temperatures, electron heat- 
ing leads to a Nernst effect, which pushes vortices away 
from the drive channel. Remarkably, this happens irre- 
spectively of the sign of the drive current in both cases. 
The qualitative features as well as the absolute values 
of the observed nonlocal voltages agree well with the re- 
sults of our model calculations. Our results offer a new 
possibility to probe the presence of vortices or vortex- 
like excitations as currently discussed in the context of 
cuprate superconductors 0. 
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